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' Abstract. We discuss results of analyzing a time series of selected photospheric-optical spectra of 

^ I core collapse supernovae (CCSNe). This is accomplished by means of the parameterized supernovae 

O ■ synthetic spectrum (SSp) code "SYNOW". Special attention is addressed to traces of hydrogen at 

' early phases, especially for the stripped-envelope SNe (i.e. SNe Ib-c). A thin low mass hydrogen 

, layer extending to very high ejection velocities above the helium shell, is found to be the most likely 

■ scenario for Type lb SNe. 
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INTRODUCTION 



^ ■ Our main goals are to identify traces of hydrogen, helium and oxygen in type Ib-c SNe, 

and to identify any systematic similarities and differences correlating with other physical 



properties in the CCSNe family. 

We have selected a sample of 45 photospheric spectra of 20 CCSNe objects. In the 
present paper we show results for a few objects. A detailed analysis for the complete 
[ sample has been presented recently in Elmhamdi et al. 2006 [10]. For the purpose 

[ of our analysis we make use of the parameterized supernova synthetic spectrum code 

"SYNOW". The main task is to identify lines and put constraints on the velocities. 
^ ! SYNOW has a number of free input fitting parameters. The most important are: 1 . (Tref), 

^ I the optical depth of "the reference line". 2. (r^^): the underlying blackbody continuum 

temperature. 3. (Vphot)'- the velocity at the photosphere, estimated from the match with 
Fe II lines. The radial dependence of the line optical depths can be chosen to be either 
exponential with an e-folding velocity "Vg" as free parameter (i.e. T exp{—v/ve)), or 
a power-law with an index "n" (i.e. T v^"). Detailed description of the code can be 
found elsewhere [14, 3, 10] 

An important aspect to understand within the "SYNOW" line formation logic is the 
"detachment" concept. When assigning a minimum ion velocity, Vm,„, greater than the 
photospheric one Vphot, the ion is said to be "detached" from the photosphere, and 
consequently has a "non-zero" optical depth only starting at v„j,„. In "SYNOW" the 
profile of a detached line has a flat-topped emission, and the absorption minimum is 
blueshifted by the detachment velocity. An undetached line has a rounded emission peak. 
A slightly detached line has a flat top but only over a small wavelength interval. These 
possibilities are illustrated in Fig. l(left), where the code run adopts Vp/,of=10000 km s^^ 
and Thb=6000 K, including only hydrogen Balmer lines with T(Ha) = 10 and having 3 
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FIGURE 1. Left: the "SYNOW" run for Balmer H I lines, using Vp/,„(=10000 km s"^ and Tfe/,=6000 K. 
The resulted profiles for different detachment velocities are illustrated. Right: The SSp fit, thin line, to the 
spectrum of SN 19901 near maximum. The region around the 6250A trough is zoomed in the window. 



different detachment velocities. 

We introduce the parameter "contrast velocity", defined as: Vcont{line) = Vmin{line) — 
Vphot- We discuss as well the behaviour of a similar parameter, defined instead as a ratio, 
i-e- = Vmin{line)/vphof 



HYDROGEN ISSUE IN CCSNE 

In the following we present results for some individual objects of the CCSNe sample. 

SN 19901 : the event is one of the well sampled and exploited, both photometrically 
and spectroscopically, objects among type Ib-c SNe[9]. Fig. l(right) compares the ob- 
served spectrum at maximum light to a synthetic spectrum (SSp) that has a velocity at 
the photosphere Vphot = 12000 km s^^ and T^^ =14000 K. The introduced lines are He 
I, Fe II, Sc II, Mg II, O I, Ca II and Ha. Absorption troughs of P-Cygni He I fines at 
5876A, 6678A and 7056A are evident, T(Hel)- 2.9 and Vm,>,(HeI)= 14000 km s^^ 
although they are not simultaneously fitted (i.e. their relative strengths). This limitation 
is faced each time we analyze and fit He I lines using the LTE approach. A more precise 
analysis requires NLTE treatment as He I lines are probably non-thermally excited by 
the decay products of ^^Ni and ^^Co [15]. Apart from He I and Ha lines, the remain- 
ing lines are undetached. The absorption minimum near 6250A is well fitted by Ha, 
with v^(„(Ha)= 16000 km s~^ and assigned a moderate optical depth of 0.6. We tested 
however other plausible alternative identifications to account for a similar feature seen 
in type lb SNe, namely Si II 6355A, Ne I 6402A and C II 6580A lines. For SN 19901, 
Vconti Ha) is only 4000 km s~\ hence undetached Ne I line is excluded because it is too 
blue to fit the observed feature. Si II 6355A is obviously also too blue. C II 6580A re- 
mains then a plausible alternative for Ha, since its rest wavelength is slightly redder 
than Ha and therefore has its contrast velocity as a free parameter. In the window of 
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FIGURE 2. Left: SSp fits of SN 19901 spectra at 12d and 21d. Right: SSp fit of SN 1991D around 3 
weeks. The top panel illustrates the Ne I possibiUty, while Ha case is shown in the lower panel. 



Fig. l(right), we show the Ha and C II fit cases. C II, with Vto>i(CII)=17000 km s"^ and 
T(Cn)~ 0.003, provides a fit as good as Ha. But an acceptable C II 6580A line would 

1 ° 

require an additional velocity of about 820 km s~ compared to Ha (~ ISA difference). 
Furthermore, the He I lines have Vto„;(He I)=2000 km s^' while Vcont{Cll)=5000 km 
s^^ This means that the C II lines would need to be formed above the helium layer, 
which would be surprising, if not unphysical, in type lb SNe. SN 19901 presents there- 
fore good evidence for an Ha feature. It is not absolutely clearcut however that Ha is 
always responsible for the absorption at about 6300A seen in type lb SNe. The presence 
of the Hj8 Balmer line would be of great support. Unfortunately the optical depth suf- 
ficient to fit the Ha trough is so small that the other Balmer features are too weak to 
be unambiguously detectable (T(Ha) is about 7 times greater than T(Hj!3)). The second 
factor that makes Hj8 barely discernable is when the contrast velocity of Ha is high as 
can be clearly seen from Fig. l(left). 

Based on our investigation, the ions that generally might be encountered in shaping 
the "6000-6500A" wavelength range are: Ha, Ne I, C II, Si II, Sc II, Ca I, He I, Fe 
II, Si I and Ba II. Here, we propose some selection criteria to decide on the presence of 
hydrogen for the 6300A feature: 

1— undetached Ne I 6402A line is rejected once its too blue to fit the 6300A trough 
or/and when the other Ne I lines clearly introduce unwanted features. Similar reasoning 
applies to the Si II 6355A line. 2— with its contrast velocity as a free parameter, C II 
6580A line could be a candidate for the 6300A trough; nevertheless it is ruled out once 
it exceeds the He I contrast velocity. 

The observed spectra at 12days and 21days are compared in Fig. 2(left) to the syn- 
thetic ones with (Vphot='^0000 km s^^; Tfe^=5500 K) and (Vphot=9500 km s^^; Tift=5400 
K), respectively. He I is still detached from the photosphere with Vcom (He I)= 3000 km 
s~^ at 12d and Vconf(He I)= 2500 km s~^ at 21days. The He I 7065A seems to increase 
in strength relative to the He I lines at 5876A and 6678A, indicating the non-thermal 
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FIGURE 3. Left: SSp fit of SN 1999ex at day 4. Conspicuous line features are indicated. Right: the 
4days IR spectrum compared with SSp that contain only Unes of He I. 



excitation effects are changing but still existent. For the "6000— 6500A" wavelength 
region, we find that Fe II, He I and Ha are sufficient to reproduce the overall shape in 
the SSp at both 12days and 21days. We have, however, tested the Ne I, Si II and C II 
possibilities as candidates for the 6250A trough. They are all ruled out by means of our 
above criteria. We therefore consider Ha to be the most likely explanation in SN 19901. 

SN 1991D : the event is a peculiar case where its difficult to decide between Ha 
and Ne I. The object has narrow features and lower velocities compared to other type 
lb events at similar phases [2, 3]. Fig. 2(right) illustrates the fit to the 21d spectrum 
{vphot = 4600 km and T^^ =7000 K=. The He I lines are evident through our 
SSp (Vc.o„f(He I)= 1400 km s^^ and T(He I)=1.8), with the Na I D feature (T(NaI)=4) 
contributing to the He I 5876A broad P-Cygni profile. The presence of Na I lines 
is consistent with the good fit around 8100A. The undetached Ne I lines possibility 
is illustrated in the top panel (T(NeI)=2), while the Ha case is tested in the bottom 
panel. Except for adding Ba II lines in order to help the fit at ~6100A, the other ion 
parameters are kept unchanged. The fit to the observed absorption near 6300A with 
Ha, VcoMf(Ha)= 7400 km s^ and T(Ha)=0.46, is slightly better compared to the Ne I 
case. However the SSp, in the bottom panel, does not account for the observed features 
near 6630A and 6840A as does Nel lines in the SSp of the top panel. Ne I remains 
hence a strong candidate in this type lb event. 

SNe 1999ex: Because of weak optical He I lines, the object was classified as an 
intermediate case between lb and Ic SNe, while the evident trough at ~6250A was 
attributed to Si II 6355A [12]. 

In Fig. 3(left) the observed 4d spectrum is compared with the best fit that has Vpkot = 
10000 km s^ and T^^ =5800 K . To check the consistency of He I identification 
we combine the IR spectrum with the optical one. In fact assigning the following 
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FIGURE 4. Left: SSp fit of SN 2000H at day 19. Right: SSp fits of SN IIP 1999em observed spectra at 
day -6(top) and day 25 (bottom). The reported phases are normalized to maximum in type Ib-c SNe. 



parameters: Vcon/(He I)= 1000 km s and T(He I)=2.35, we obtain a good match to 
the observed He I profiles including the IR He I 10830A. Furthermore, unambiguous 
evidence for the presence of helium in this event is clear from the extended fit to other 
He I-IR lines, adopting equal v,„,„(He I) and T(He I) as in the optical part of the spectrum 
(Fig. 3-right). The 01 7773A line with T=0.5, on the other hand, is found to be not so 
strong and deep as is the case in most type Ic SNe [16]. These two facts point to a lb 
nature, more than a Ic class. 

As far as the trough around 6250A is concerned, we checked the Si II identification 
attributed by Hamuy et al.(2002)[12]. We test as well the Ne I possibility. The unde- 
tached Ne I fit introduces unwanted features in the spectrum. The Si n 6355A is rather 
blue to account for the feature. We note here that if one adopts C II 6580A, then one 
needs to assign it a very high velocity, about 8000 km s~^ higher than the one of He 
I. The most likely identification therefore remains Ha. In fact the best fit is achieved 
using Vcon?(Ha)= 8000 km s~^ and T(Ha)=1.45, 

SN 2000H : the SN is an example of a type lb event with a strong and deep 
trough near 6300A [3]. Fig. 4(left) compares the 19d spectrum with the best fit SSp 
{vpi^„f = 6000 km s-i and Tbb =4600 K). Lines of Ca H, [O II], O I, Mg H are seen 
to develop. The optical He I lines are noticeably deep (vcoMr(He I)= 2000 km s^^ and 
T(He I)=5). We note the good match to the ~6250A trough with Ha, achieved using 
Vconr(Ha)= 7000 km s^^ and T(Ha)=2.5. No alternative identification to Ha, that would 
be logically acceptable, has been found. In addition, the Ha identification is supported 
by the presence of the absorption notch near 4660A, well accounted for by Hj8 in our 
resulting synthetic spectra. 

SN IIP 1999em: Fig. 4(right) displays two spectra during the photospheric phase 
of the typical type Il-Plateau SN 1999em [8]. The phases are normalized to maximum 
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light in type Ib-c SNe. The best fit synthetic spectra have Vphot = 10000 km s and 
T^,^ =10000 K (-6d; upper panel) and Vphot = 4600 km s'^ and T^^ =6400 K (25d; 
lower panel). At the earliest phase only undetached lines of hydrogen B aimer. He I 
and a weak contribution from Ca n are sufficient to reproduce the most conspicuous 
features superimposed on the "hot" continuum. The fit with Na I D, at this phase, is poor 
compared to He I 5876A. The He I contribution, in type II SNe, is found to be important 
shortly after the explosion. In fact in SN 1987A, the He I 5876A feature was clearly 
present during the first few days, then rapidly faded and disappeared completely around 
1 week after the explosion, when Na I D starts to emerge [13]. At -6d, Balmer hydrogen 
features are evident with T > 15 for Ha. The corresponding strong and broad P-Cygni 
feature cannot be produced completely by the "SYNOW" code. This is because in our 
SSp treatment we are adopting a resonance scattering source function. At day 25, Fig. 
4(right-bottom), Ha is narrower and has Vcont= 1000 km s~^ and T = 21. The envelope 
temperature decreases and many lines emerge at this phase. Apart from the hydrogen 
lines (slightly detached), all the lines are undetached. 



RESULTS AND DISCUSSION 

Based on results presented here and in the larger sample study of Elmhamdi et al. 
(2006) [10], hydrogen is found to manifest its presence in a different way in the CCSNe 
sample. SNe of type n and lib show strong and broad Ha P-Cygni emission component, 
which is absent in type lb objects. Indeed, the Ha P-Cygni profile would lose its obvious 
emission component when it is highly detached. H I Balmer lines have large optical 
depths, allowing them to be distinctly visible, contrary to normal type lb SNe, with 
the exception of some cases with deep and conspicuous Ha troughs that present the 
signature of H/3 as well. Two factors make H/3 barely discernible in type lb: the optical 
depth found to fit Ha is small and the contrast velocity of Ha is high. The opposite is 
seen in type II and lib events. 

In Fig. 5, upper panel, we report the resulting photospheric velocities from our best 
fits for the whole CCSNe sample [10]. Data for SN 1987A, corresponding to Fe n 
5018A absorption, are also shown for comparison (dashed line [17]). The plot indicates 
the low velocity behaviour of type II SNe, both 1987A and 1999em. SN lib 1993 J 
follows somewhat similar behaviour as SN Ic 19941 in having higher velocities. Type 
lb SNe appear to display a different velocity evolution. The scatter seems to increase at 
intermediate phases (around 20—30 days). This fact can be simply due to the paucity 
of available observations outside that range. Around day 20, for example, a scatter as 
high as 5000 km s~^ is measured. SNe 19901 and 1998dt belong to a class with the 
higher ''vpho", while objects such as SNe 1991D and 1996aq have the lowest estimated 
velocities, approaching even type II objects. The remainder of the lb events follow a 
similar trend, namely the one described by Branch et al. (2002) [3]. 

The middle and bottom panels in Fig. 5 display the evolution the of Ha contrast 
velocities. An increasing trend of "vconf(Ha)" is found going from type II to lib to lb 
SNe. While in SN lib 1993 J and SN IIP 1999em the line is found to be either undetached 
or slightly detached, it is highly detached in type lb events. Moreover, the "vcon? (Ha)" 
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FIGURE 5. The photospheric velocity evolution of the analyzed CCSNe (top). The Ha contrast velocity 
evolution is also illustrated (middle and lower panels). 



is found to increase within the first 15 days, reaching values as high as 8000 km s~ , 
and then follows an almost constant evolution. Up to ~60 days, type lb SNe are found 
to have hydrogen down to 1 1000-12000 km s~^ SNe n appear to have hydrogen down 
to lower velocities (~5000 km s^^ in SN 1999em). In addition, hydrogen in type lb SNe 
is found to have very small optical depths independently of the contrast velocity. 

An approximate method to estimate the hydrogen mass from the fit results invokes 
the amount needed to fill a uniform density sphere of radius "v x f at an epoch t^j since 
explosion. The required ion mass can be given by ([10]): 

M(M0)~ (2.38 X 10"^) v|rjT(//a); (1) 

where in days and V4 is in 10"^ km s ^ Although non-thermal excitation and "NLTE" 
effects may be also important for hydrogen, this simple approach seems to give rea- 
sonable estimates. For SN lb 19901 around maximum light a value of 0.02 Mq is com- 
puted. If we adopt a representative Ha optical depth of 0.5 at day 20 with an Ha ve- 
locity restriction similar to SN 19901, the estimated hydrogen mass is of the order 0.015 
Mq. Events with higher velocity widths and/or deeper Ha troughs would eject larger 
amounts. A thin layer of hydrogen, ejected at high velocities down to 1 1000-12000 km 
8"^ appears to be present in almost all the studied type lb events. 

Worth-noticing here is that the presence of Ha in type Ic events was recently re- 
addressed by Branch et al. (2006) [4]. A large type Ic sample study is clearly needed 
to assess this issue. However, based on these preliminary investigations the progenitor 
nature of the stripped-envelope SNe should be revised. 

An important by-product results of our fit analysis concerns the behaviour of the O I 
7773A line. At intermediate phases, it seems that type lb objects tend to have low optical 
depths in this line, while SNe Ic are found to display the deepest profile. SNe lib an II, at 
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similar phases, are the objects with the lowest 1 1113 A optical depth. The stronger and 
deeper permitted oxygen lines at early phases of SNe Ic spectra might imply that they 
are less diluted by the presence of a helium envelope. Indeed one might expect oxygen 
lines to be more prominent for a "naked" C/O progenitor core. Despite the paucity 
of well sampled CCSNe observations, two observational aspects tend to reinforce this 
belief: First, the forbidden lines, especially [OI]6300, 6364A, seem to appear earlier 
following a SNe sequence "Ic— lb— lib— 11". In fact the oxygen line emerges at an age of 
1-2 months in Type Ic SN 1987M [11]. SN Ic 19941 displayed evidence for the line at 
an age of 35-50 days [7]. While in SN lb 19901 it was hinted at the VOday spectrum 
[9]. In other Type lb SNe it appears earlier than in SN 19901. In SN lib 1993J, a 
transition object, the line was visible in the 62day spectrum [1]. In SNe II, however, 
the line appears later: around day 150 in SN 1987A [5] and after day 138 in SN 1992H 
[6]. Second, the nebular emission line looks decreasing in breadth following the SNe 
sequence above. More investigations concerning the oxygen issue in CCSNe are under 
progress (Elmhamdi & Danziger, in preparation). 

REFERENCES 

1. Barbon R., Benetti S., CappeUaro E. et al. 1995, A&A5, 110, 513 

2. Benetti S., Branch D., Turatto M. et al. 2002, MNRAS, 336, 91 

3. Branch D., Benetti S., Kasen D. et al. 2002, ApJ, 566, 1005 

4. Branch D., Jeffery D.J., Young T.R. et al. 2006, " astro -ph/0604047" 

5. Catchpole R. M., Whitelock R A., Feast M. W. et al. 1988, MNRAS, 231, 75 

6. Clocchiatti A., Benetti S., Wheeler J. C. al. 1996a, AJ, 111, 1286 

7. Clocchiatti A., Wheeler J. C, Brotherton M. S. et al. 1996b, ApJ, 462, 462 

8. Elmhamdi A., Danziger 1. J., Chugai N. N. et al. 2003, MNRAS, 338, 939 

9. Eknhamdi A., Danziger I. J., Cappellaro E. et al. 2004, A&A, 426, 963 

10. Elmhamdi A., Danziger I. J., Branch D. et al. 2006 , A&A,450 , 305-330 

11. Filippenko A. V. 1997, A/M&A, 35, 309 

12. Hamuy M., Maza J., Pinto R A. et al. 2002, AJ, 124, 417 

13. Hanuschik R. W. & Dachs J. 1988, A&A, 205, 135 

14. Jeffery D. J. & Branch D. 1990, in: J.C. Wheeler, T. Piran, & S. Weinberg (eds.),"Supernovae", 
Sixth Jerusalem Winte School forTheoretical Physics (Singapore: World Scientific), 149 

15. LucyL.B. 1991, Ap/, 383, 308 

16. Matheson T., Filippenko A. V., Li W., Leonard D. C. & Shields J. C. 2001, AJ, 121, 1648 

17. PhilUps M. M., Heathcote S. R., Hamuy M. & Navarrete M. 1988, A7, 95, 1087 



Hydrogen issue in Core Collapse Supemovae 



February 5, 2008 



8 



